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Abstract 

This work proposes a new way for handling obstacles to asymp- 
totic integrability in perturbed nonlinear PDE's within the method 
of Normal Forms (NF) for the case of multi-wave solutions. Instead 
of including the whole obstacle in the NF, only its resonant part 
(if one exists) is included in the NF, and the remainder is assigned 
to the homological equation. This leaves the NF integrablc and its 
solutions retain the character of the solutions of the unperturbed 
equation. 

We exploit the freedom in the expansion to construct canonical 
obstacles which are confined to the interaction region of the waves. 
For soliton solutions (e.g., in the KdV equation), the interaction 
region is a finite domain around the origin; the canonical obsta- 
cles then do not generate secular terms in the homological equation. 
When the interaction region is infinite (or semi-infinite, e.g. in wave- 
front solutions of the Burgers equation), the obstacles may contain 
resonant terms. 

The obstacles generate waves of a new type which cannot be 
written as functionals of the solutions of the NF. When the obstacle 
contributes a resonant term to the NF, this leads to a non-standard 
update of the wave velocity. 

Keywords; nonlinear evolution equations, wave interactions, obstacles 
to asymptotic integrability, perturbed KdV equation, perturbed Burgers equa- 
tion. 

1 The problem of obstacles to asymptotic inte- 
grability 

The analysis of the effect of a perturbation on wave solutions of evolution PDE's 
has evolved in two different approaches. In scattering approach, a solution of 
the unperturbed equation is scattered off a perturbation that is turned on at 
t = 0. As the unperturbed solution is not a solution of the perturbed equation, 
its amplitude decays and its wave number, velocity and phase shift are modified. 
In addition, a soliton tail has been found outside the soliton sector in the case 
of the perturbed KdV equation. The methods used have been a combination 
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of the theory of inverse scattering and the multiple time expansion procedure 

The second approach deals with a different issue. Integrable nonlinear evolu- 
tion equations are lowest-order approximations for the more complicated equa- 
tions of the full dynamical systems (e.g., the equations of Fluid Dynamics in 
the cases of the Burgers & KdV equations, and Maxwell's equations in the case 
of the NLS equation) . To improve the approximation, one includes higher-order 
effects of the original physical system. In this case the perturbation is not turned 
on at t = 0. It exists for all times, going back to t = — oo. 

In this approach, one is searching for a zero-order solution, which has the 
same wave structure as the solution of the unperturbed equation, except for 
an update of the wave velocity by the higher-order effects. The method of 
Normal Forms was used for the analysis of solition solutions of the perturbed 
KdV and NLS equations H2]-[I3|, CD- [22] and front solutions of the perturbed 
Burgers equation JH], 123- The me thod of Multiple Time Scales was used 
in the case of the perturbed KdV equation ^8], |2J and in the case of the 
perturbed NLS equation |2.'{| . 

We focus on problems that arise in the analysis of wave solutions within the 
second approach. Most integrable nonlinear evolution PDE's allow for single- 
as well as multi-wave solutions. The multi-wave solutions usually asymptote to 
well-separated single waves in the x — t plane, except for interaction regions, 
where the multi-wave character of the solution is lost. The interaction regions 
may be localized (e.g., in the case of KdV-multi-soliton solutions) or semi-infinite 
(e.g., Burgers-multi-fronts). The main purpose of this work is to investigate the 
effect of a perturbation that is added to the unperturbed equation on the wave 
solutions of nonlinear systems. 

The perturbed equations are often analyzed by the method of Normal Forms 
(NF) [5] - |TI], briefly described in the following. Let 

w t =F(°'H+^e fc Ft l )[io] (1) 
fc=i 

be a perturbed nonlinear evolution PDE (square brackets imply that the 
corresponding term is a differential polynomial in w(x, t)). We assume that 
w (x,t) may be expanded in a power series in the small parameter e of differential 
polynomials of u (x, t) (NIT - Near- Identity Transformation): 

w=J2^ u{k) [u] (u = u (0) ) (2) 

fe=0 

The time evolution of the zero-order term, u{x, t), is assumed to be governed 
by the Normal Form (NF): 

u t =JVa fc SW[u], (a = l) (3) 

k=0 



2 



Here, are the resonant terms usually called the symmetries. Their time 
dynamics is equivalent to that of u (x, t) up to first order, that is, 




As a result, their Lie Brackets vanish: 




The symmetries (including itself) form hierarchies |3l)|. |31|. It is possi- 
ble to establish recursion relations among the symmetries in each hierarchy. For 
many equations (and all the equations our work deal with), the first symmetry 
is 



Substituting the NIT © and the NF © into (JTJ leads to a sequence of 
homological equations for the time evolution of u^ n \ which have to be solved 
order-by-order . 

The motivation for assumption (3) is that a perturbative analysis that does 
not include the removal of resonant terms from the homological equation into 
the NF, usually yields secularities, that is, unbounded terms in the approximate 
solution. On the other hand, the NF is expected also to be integrable and to 
preserve the nature of the unperturbed solution. This feature is closely related to 
another significant one: that the main effect of adding the higher-order terms 
to the NF is the update of physically-valuable parameters (usually the wave 
velocity /the dispersion relation). 

After removing the resonant terms out of the homological equations, they 
become: 



where \ u] is the contribution for all nonresonant terms of order k. The NIT 
is constructed from solutions of these equations. 

However, the analysis may lead to the emergence of obstacles to integrability 
22]- E3 j beginning at some order in the expansion. These are terms (differential 
polynomials) that the perturbative expansion of the dynamic equation (Eq. ^| 
generates, which cannot be accounted for by the formalism. The differential 
polynomial structure of the obstacles is not unique and depends on the way in 
which the NIT is constructed. 

To make the construction of the NIT possible, the usual practice has been 
to include these unaccounted-for terms in the NF. This makes the NF nonintc- 
grable, hence the name "obstacles to integrability" . Including the obstacles in 
the NF, disturbs the wave character of its solutions. The effect of the obsta- 
cles in the case of the two-soliton solution of the normal form of the perturbed 
KdV equation 22-22 nas been studied in 122 . The zero-order solution was 




(G) 



u (k) [u] + T< k > [u] = 



(7) 
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found to develop inelastic effects: appearance of a second-order radiation wave; 
fourth-order, time dependent, corrections to each of the wave numbers and the 
generation of an eighth-order soliton. 



2 Approach for overcoming obstacles 
2.1 The general ideas 

The necessity to include the obstacles in the normal form is a consequence of 
the assumption, usually made in the NF expansion, that all the terms in the 
NIT are differential polynomials in the zero-order approximation (that is, in 
u(x,t)) and do not depend explicitly on the independent variables, t and x. 
Our approach overcomes this problem by allowing the higher-order terms in the 
NIT to depend on these variables. To this end, we assume for the fc'th order 
term in the NIT the following form: 



U W = „<*> [u] + 



(x,t) (8) 

In Eq. |JBJ, u'p [u] is a differential polynomial in u, and (x,t) depends 
explicitly on x and t, and is expected to account for the obstacles. Thus, sub- 
stituting the assumption (JHJ in the homological equation JJJ , we obtain: 



F ( ' , u ( ;>(x,t) +R {k >[u}=0 (9) 

where [u] stands for the obstacle of order k. 

Owing to the freedom inherent in the perturbative expansion, the construc- 
tion of [u] is not unique. Unless [u] is chosen in an appropriate manner, 
the resulting obstacle may not reflect the following features of physical interest: 

(i) Obstacles do not emerge in the case of single-wave solutions of the NF 

m 

(ii) The expectation that obstacles emerge owing to interaction among waves 
in the multi-wave case |17|. 

Both features are realized if Ug [u] is chosen to have the structure of the 
differential polynomial that solves the problem in the case of a single-wave 
solution of the NF. The choice proposed for uf^ [u] leads to obstacles in a 
"canonical" form, expressed in terms of symmetries of the unperturbed equation. 
The obstacles now vanish if one substitutes for u the single-wave solution of the 
NF. More important, as a result, they are expected to vanish away from regions 
of wave interaction in the multi-wave case. The reason is that away from the 
interaction regions, multi-wave solutions asymptote into a sum of well-separated 
single- wave solutions. 
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2.2 Construction of canonical obstacles 

Starting with [u] that solves the homological equation J7J) for the single- 
wave, we find that the canonical obstacles can be written in the following form 

Rto [ u ] = 7 fc fij [u, d x ] Rij [u] (10) 

fe=3 
i+j=k 

R pq [u] = S« [u] G« [u] - S<«> [u] GW [u] , 

£(«) [ u ] = d x G (n) [u] [ ' 

where 7J? is a numerical coefficient, fij [u, d x ] - a differential operator of an 
appropriate weight and g is the gap between the index of symmetry and the 
order of perturbation. For example, if we define as (this is a widely- 
accepted notation), then g = 2. 

The obstacles of Eq. I|10|) vanish identically for the case of single-wave so- 
lutions of the NF. To see this, we exploit the fact that all the symmetries are 
proportional to one another in the case of the single-wave solutions of the NF. 
For "trivial" boundary condition, u (£ — ► -co) = 0: 

S< n > =(-l) n+1 (12) 

This may be proven by induction for any hierarchy which is governed by a linear 
recursion relation. 

The proportionality of all the symmetries leads to a simple update of the 
velocity of the solutions of the NF: 

v = ^ ■ v k = ("l) fe a k v k +1 (13) 

fc>0 

Eq. (|13|) also describes the velocity update of each wave in a multi-wave solution. 

2.3 Resonant contribution in obstacles? 

For multi-wave solutions of the NF, the obstacles do not vanish. An important 
characteristic of our canonical obstacles is that they do not vanish only in the 
interaction regions in the x — t plane. For example, in the case of KdV solitons, 
the interaction region is a finite domain around the origin, whereas in the case 
of Burgers fronts it consists of one ore more domains of finite width along 
semi-infinite lines. The canonical obstacles vanish exponentially fast away from 
the interaction regions, where the solution asymptotes to well-separated single 
waves. On the other hand, non-canonical obstacles are finite also outside the 
interaction region. 

A cardinal question that now arises is whether the obstacles generate secular 
terms in the NIT (J2J. A symmetry, if contained in an obstacle, will generate 
a secular term through Eq. 10. We, therefore, propose to break an obstacle 
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into a sum of a symmetry plus a non-resonant term. This symmetry, with its 
coefficient, must be included in the NF ©. 

Our task is therefore to determine whether a canonical obstacle has the 
capacity of generating secular terms. A simple criterion for detecting this ca- 
pacity is that the obstacle spreads over an infinite or a semi-infinite domain, 
and asymptotes to a symmetry. (This criterion is similar, although much less 
rigorous, than the Fredholm Alternative Theorem.) 

We focus on two-wave solutions of two equations: the Korteweg-de Vries 
(KdV) equation with a two-soliton solution, with obstacle appearing in the 
second order; and the Burgers equation, with a two- front solution, where an 
obstacle arises already in the first order. 

The interaction region, and thus the obstacles, of KdV are localized. Hence, 
we expect the solution of Eq. [§]for this problem to be bounded. This expectation 
has been verified by solving the homological equation numerically. We obtain a 
new bounded soliton-like wave. 

On the other hand, the interaction region in the two-front solution of the 
Burgers equation is semi-infinite (the fronts are well separated in one half of the 
x — t plane and merge (interacting) in the other half). The canonical obstacle 
asymptotically approaches the symmetry [u]. Therefore, we expect the 
Burgers obstacle to generate a secular term in Eq ©. This expectation has 
been verified by numerically solving the equation. 

We "extract" the symmetry [u] out of the obstacle with its coefficient, 
and transfer it into the NF J2J. Thus, the NF remains solvable, but the co- 
efficient of the first-order term in the wave-velocity update is changed. The 
remainder is not a canonical obstacle, that is, it is not confined to the inter- 
action region but rather spreads over all the fronts of the zero-order solution. 
However, it does not contain a symmetry. The solution of the homological 
equation must be bounded. There is no closed form solution of the homological 
equation, but the numerical calculations show that this prediction is verified. 

In the following section, these two examples are brought in some detail. The 
full results of our work will be brought in further publications. We stress that 
our results crusually depend on the wave nature of the solutions. 

2.4 Summary and Conclusions 

1 . The main effect of perturbations on the interaction among waves in multi- 
wave solutions of NL EPDE's is the emergence of the obstacles to integra- 
bility. 

2. The proper place for handling obstacles, once their resonant part has been 
shifted to the NF, is the homological equation. To this end, it is necessary 
to allow the higher-order terms in the NIT (0 to be explicit functions 
of the independent variables. Actually, will consist of two parts: the 
differential polynomial, [u] , whose structure corresponds to the case of 
a single- wave solution of the NF © (in this case, there are no obstacles); 
and the function u£ t) that is supposed to account for the obstacle. 
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3. When the NIT is built in such a way, the obstacles (if they exist) are 
obtained in the canonical form l|10f) . Canonical obstacles are confined to 
the region of interaction among waves. 

4. When the interaction region is localized (e.g., for KdV solitons), the canon- 
ical obstacles do not generate secular solutions of the homological equa- 
tions The solution of the homological equation is then bounded, 
usually unavailable in closed form. 

5. When the interaction region spreads over a semi-infinite range, the ob- 
stacles are expected to cause secular solutions. Then, it is necessary to 
identify a symmetry "hidden" inside the obstacle and remove it from the 
homological equation into the NF. The remainder of the obstacle remains 
in the homological equation and yields bounded solutions. 

In both situations, the obstacles cease to be obstacles to integrability of the 
NF. The latter remains integrable, and its solutions (the zero-order approxima- 
tion) retain the character of the unperturbed solutions. The difference between 
the items 4 and 5 above is that in # 5 the wave velocity is also affected, in the 
order in which the obstacle exists. Focusing on two-wave solutions, it is found 
that the obstacle usually yields an additional wave, that may not be expressed in 
closed form. In general, the homological equation has to be solved numerically. 



3 Two worked examples 

3.1 The perturbed KdV equation 

The perturbed KdV equation is 

wt — 6ww x + + e (30aiw 2 w x + I0a2ww xx + 20a3W xW xx — CC^W§ 
2 / U0Piw 3 w x + 7Q(3 2 w 2 w xxx + 280{3 3 ww x w xx + 14/3 4 ww 5: 
I +70f3 5 w x + 42f3 6 w x w ix + 70(3 7 w 

XX^XXX 

We assume the NIT 




w = u + eu {1) +eV 2) +0 (e 3 ) (15) 

and the NF 

u t = S (2) [u] + ea 4 S (3) [u] + e 2 /3 8 S (4) [u] + O (e 3 ) (16) 

where 



S" (2) [u] =6uu x + u xxx 
S(3) [ u ] = 30u 2 u x + lOu 

u xx ~\~ 20u x U xx -\- U^ x 

[u] = 140u 3 u x + 70u 2 u xx + 280u x u xx + I4uu 5x 
+70u 3 + A2u x u 4x + 70u xx u xxx + u 7x 



(17) 
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A single-wave solution of the NF (| 1 fcil) is the well-known KdV soliton 

2k 2 

cosh 2 [k (x — vt + xq)] 
and the two- wave solution is given by the Hirota formula |26) : 



u 0> *) = . ~ , -i (18) 



u(x,i) =252 In{l + 5 i(a:,t)+fl2(i > t) + (*J=^) 5 i(a; > t)92(a:,t)} 
(ffi (x, t) = exp [2fcj (x - + x ,i)]) 



(19) 



A sample of this solution is shown in Fig. 1. The only difference between the 
solutions of the unperturbed KdV equation, wt — 6ww x + w xxx , and the both 
solutions of the NF is the update of the wave velocity: 

Vi = -4fcf - 16ea 4 fci - 64e 2 /3 8 fc 4 6 - O (e 3 ) (20) 

There are no obstacles in the first order in the Normal Form analysis. How- 
ever, an obstacle appears at the second order. Choosing the second-order term 
in the NIT to be 

u ( V = uf' (x, t) + Bxu 3 + B 2 u u xx + B z u 2 x + B A u xxxx 
+B 5 u u x q^ + B 6 v, xxx Q W + B 7 u xx q^ 2 + BgU x qW ( 21 ) 

with the appropriate set of values for {Bk} based on the form of w 2 ' in the case 
of the single-soliton solution, we obtain the canonical obstacle 

R = = 73W (3u 2 u x + uu xxx - u x u xx ) (22) 

(73 is built of a combination of coefficients of Eq. (|14J) . This obstacle is localized 
(see Fig. 2) and hence cannot generate a secular solution in the homological 
equation 

8 t u {2) = 68 x (u uty + dlu^ + 1 luR 2l (23) 

This equation can be solved in closed form by the Green's function method 
developed in the context of the Inverse Scattering approach |2J, |32|-|34|, As 
our goal is only to show that the solution of Eq. I|23|l , a numerical solution was 
sufficient for our purpose. It shows a new bounded soliton- like wave (Fig. 3). 

3.2 The perturbed Burgers equation 

The perturbed Burgers equation is given by 



w t = 2ww x + w xx + e (Zaiw 2 w x + 3a 2 ww xx + 3a 3 w 2 x + a±w xxx ) + O (e 2 ) (24) 
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(we don't write here the second-order perturbation because an obstacle emerges 
already in the first order). We assume, again, the NIT 



w = u + e«W + O (e 2 ) (25) 

and the NF 

u t — 2u u x + u xx + e/z (3u 2 u x + 3u u xx + 3u 2 x + u xxx ) + O (e 2 ) (26) 
Its single- wave solution is the shock front 

_ kAexp [k (x-vt)} 
U(X ' >~ l + kAexp[k(x-vt)} 1 ' 

and the two- front solution is a straightforward extention of Eq. (|27H : 

k\A\ exp [fci (x - vit)} + k 2 A 2 exp [k 2 (x - v 2 t)} fno . 

U(X,t) = ; r- ; tt - ; — (28) 

1 + hAx exp [fci (x - wit)] + k 2 A 2 exp [k 2 (x - v 2 t)] 
(its sample is shown in Fig. 4, and one can see that the fronts are interacting 

(i.e. merged) over a semi-infinite region). The velocity update is, again, similar 
for the both solutions: 

v t = -h - e^k 2 - O (e 2 ) (29) 

The obstacle appears in the first-order analysis. If we move all the linear 
term u xxx from the first order into the NF, that is, take fj, = a^, and also choose 
u' 1 ) in the form that solves the homological equation for the single-front case 
plus an explicit function of x and t: 



uM = Ur (x, t) + [ax — 2a 2 — as + 2014) qu x — | (2a\ — a 2 + — 2a 4 ) u 2 

(q = d^u) 

(30) 

then the obstacle will be get its canonical form: 



m = liiia! = 73 1 (S^G^l - SWgW) = 71 (uS< 2 ) - u x G^) 
(73 = 2a\ —012 — 2a 3 + a 4 ) 



(31) 



This obstacle is shown in Fig. 5, and one can see that it is finite over all the 
interaction region. The homologicals equation now reads 

cK (1) = 2d x (u u^) + d 2 u^ + jlR 21 (32) 

It should be remarked that of the two terms of R 2 \, only uS^ may not 
be accounted for by the differential polyomials in the NIT. On the contrary, it 
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is easy to see that substituting Ur — into the homogeneous part of the 
homological equation yields 

-d t G {2) + 2d x (u G (2) ) + d 2 x G {2) = 2u x G [2) (33) 

As expected, the numerical solution of the homological equation shows the 
existence of a secular term, which indicates the presence of a symmetry inside 

-^21- 

In order to extract this symmetry, we use the recursion relation among sym- 
metries in the Burgers hierarchy: 

S (n+1) = gin) + u gin) + ^(n) = gjp) + ^ + 2u , x G^ (34) 

In particular, for n = 2, it reads 

5 (3) = S {2) + R 2 i + 2u x G (2) => R 2 i = S {3) - S x 2) - 2u x G {2) (35) 
and the homological equation l|32l) becomes 

d t u^ = 2d x + flguW + 73 1 (s< 3 > - S X V - 2u x G^) (36) 

Now, we choose fi, the coefficient of in the NF, to be 

fi = Q4 + 73 = 2ot\ — a-i — 2a3 + 2«4 (37) 

hence correcting the update of the wave velocity (Eq. I29|) . Further, we add an 
appropriate correction to in order to account for the term — 2"f^u x G^ in 
Eq. (5Bjl. according to Eq. (p^ . 

Now, the equation for ui- 1 ' becomes 

8 t u^ = 2d x (u u«) + d 2 x u^ - 73^ (38) 

(2) 

The term S x does not asymptote to any symmetry. The numerical solution 
of Eq. H38f) shows that a new bounded wave appears (see Fig. 6). 
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5 Figures 

Figure captions: 

Fig. 1: Two-soliton solution of the KdV NF (Eq. UU; ki = 0.3, k 2 = 0.4. 
Fig. 2: Canonical obstacle wi?2i for the KdV equation fEa l22f) for two- 
soliton solution; parameters as in Fig. 1. 
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,(2) 



of Eq. (EU for 



Fig. 3: Contribution of this canonical obstacle to uf 
zero boundary condition for x — > — oo; fei = 0.5, k% = 0.75. 

Fig. 4: Two-front solution of the Burgers NF (Eq. EBJ; fci = 2, fc 2 = -2. 

Fig. 5: Canonical obstacle — i?2i (with opposite sign) for the Burgers 
equation (Eq. I31fl for the two-front solution; parameters as in Fig. 4. 

Fig. 6: Bounded contribution of the obstacle to of Eq. (|3*5|) : 
parameters as in Fig. 4. 
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Fig. 1 



Fig. 2 
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Fig. 5 



Fig. 6 
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